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SUMMARY 


Owing to the bulk of the material presented, this 
report is divided into several volumes. Volume 1 introduces 
the series and presents and discusses the results for Shot 10. 
Subsequent volumes will present and discuss the results for 
Shots 8, 9 and 11. The method of analysis is common to all 
four experiments and is described in detail only in Volume l. 

So that the results from the four experiments may be 
easily compared, they have been scaled to remove the effects 
of varying atmospheric conditions. (Results are scaled to a 
1 kg charge weight and a standard atmosphere of dry air at 
15 °C at sea level.) For the most part, only scaled results 
are presented. Exceptions include some derived pressure-time 
histories, which are meant to be compared to actual gauge 
measurements made in each of the experiments. 

Results are presented in SI units, even though the 
experiments were originally laid out in British units. Only 
distance and time measurements are affected, however, as velocity, 
density, and pressure results are presented as dimensionless 
ratios. A distance units conversion scale is included for 
convenience, to convert between SI units (meters scaled to a 
1 kg charge) and British units (feet scaled to a 1 1b charge), 


plus a time scale factor and scale factors to convert pressure 


ratios to both British and SI pressure units. Scale factors 


which may be used to compute the distance and time values 


actually observed under the ambient conditions of each 
shot are also provided. Real pressure units are used for 


the results presented at gauge locations. 
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CHAPTER 1, INTRODUCTION 


1.1 Blast wave reflection 
When a spherical charge is detonated above the ground 
surface the resulting spherical blast wave reflects from 
the ground. At a distance from the point on the ground 


immediately beneath the charge, approximately equal to the 


height of the charge above the ground, the reflected shock 
begins to overtake and combine with the primary shock to form 
a Single shock known as the Mach stem. The point at which 
the primary shock, the reflected shock and the Mach stem meet 
is called the triple point. As time progresses and the Mach 
stem shock moves outwards, the triple point moves outwards and 
upwards in a curved trajectory. 

The physical properties of the Mach stem blast wave and 
the trajectory of the triple point depend primarily on three 
things: the energy yield of the charge, the height of the charge 
above the ground, and the nature of the ground surface. When 
the primary shock reflects from the ground some energy will be 
absorbed and appear as seismic disturbances, including cratering 
if the explosion is close enough to the ground. As the blast 
wave continues to move across the ground surface there will be 


a continued transfer of energy between the air and the ground, 


receding Fege BLavh - 
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and also a redistribution of energy within the blast wave. 
Little is known about the transfer of energy to the ground 
from an air blast, or about the redistribution of energy 
within a blast wave as it passes over the ground surface. 

In the experiments described in this report, attempts 
have been made to simulate the "ideal" reflection of spherical 
blast waves to provide a reference for studies of reflection 
by real ground surfaces. It is postulated that if two 
identical charges are simultaneously detonated at a certain 
distance apart, the two resulting identical spherical blast 
waves will interact along a plane, and since there will be no 
energy loss in this interaction it will be possible to observe 
an ideal spherical blast reflection. The properties of the 
Mach stems lying above and below the ideal reflecting plane 
may then be compared with those produced over various real 


ground surfaces. 


1.2 General description of the project 


The primary purpose of the project was to obtain informa- 
tion on the interaction of spherical blast waves from explosive 
sources with real and ideal reflecting surfaces. The results 
may be used, for example, to evaluate hydrodynamic air blast | 
computer codes. The blast wave interactions were obtained by 


the simultaneous detonation of two identical spherical charges 


he 


placed one above the other, such that the distance between 
the charges was equal to twice the height of the lower charge 
above the ground surface. 

Two different charge heights were used over two different 
types of ground surface. The four experiments were titled 
Dipole West Shots 8, 9, 10, and 11. For Shots 8 and 11 the 
lower charge was placed at a height of 25 ft above the ground 
surface and the second charge at an additional height of 50 
ft above the first. Shot 8 was carried out over smooth ground 
and Shot 1l over rough ground. For Shots 9 and 10 the correspond- 
ing figures for charge height and charge separation were 15 ft 
and 30 ft. The ground surface was smooth for Shot 9 and rough 
for Shot 10. 

Two photogrammetrical studies were carried out in each 
of these experiments. The first study involved the high 
speed photography. against a suitable background, of the shock 
fronts produced by the two explosions. This permitted a 
calculation of the shock front trajectories and velocities, 
and thus a determination of the physical properties immediately 
behind the shocks as they moved into ambient air. The results 
of those calculations have been reported by Dewey et al (1975). 
The second project involved the high speed photography of an 
array of smoke puffs which acted as particle flow tracers, 


to determine the particle trajectories within the blast waves. 
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These trajectories have been analyzed to provide the space and 


time variations of particle velocity, density, hydrostatic 
overpressure, and dynamic pressure within the waves. 

This volume describes the photogrammetry of the particle 
trajectories for Shot 10. The particle trajectory photo- 
grammetry for the other three experiments will be reported 
in subsequent volumes. Shot 10 was analyzed first because 
the data collected were fewer than those collected for 
Shots 8 and 11 (fewer smoke puffs were used in Shot 10), and 
of better quality than those collected for Shot 9 (all smoke 
puffs detonated in Shot 10 and most remained visible throughout 
the time interval of interest). Thus the analysis of Shot 
10 was carried out with a minimum of data processing and 


maximum expectation of useful results. 


a 


1.3 Description of Shot 10 
Dipole West Shot 10 was fired on November 2nd, 1973 by 


the Ballistics Research Laboratories at the Defence Research 
Establishment Suffield, in Alberta, Canada. Two 1080-lb. 
spheres of Pentolite were detonated simultaneously, to within 
5 microseconds, at nominal charge heights of 15 and 45 feet 
over a relatively rough (ploughed) ground surface. 

Particle trajectory data were gathered by photographing 


the movement of smoke puffs formed in a vertical plane running 
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out from ground zero at 6.7° south of west. A WF5 camera 
Operating at about 3400 frames per second was positioned 30 
ft above the ground level at the main camera position 600 ft 
due south of ground zero. 

Figure 1 shows a plan view of the field layout. The 
dashed line represents the approximate line of sight of the 
WF5 camera. Figure 2 shows the field of view of this camera. 

The smoke puff grid was made up of 9 columns of 12 
puffs each, hung vertically on strings. The vertical spacing 
of puffs was 5 ft, beginning 3 ft above ground level and ending 
at a height of 58 ft. The horizontal spacing of the columns 
of puffs was 10, 7 or 5 ft, depending on the distance from 
ground zero, beginning at about 25 ft and ending at about 85 
ft. All 108 smoke puffs detonated successfully and a good 
film record was obtained. 

This report describes the analysis of the smoke puff 
data collected for Shot 10, and presents and discusses some 


of the results of that analysis. 
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CHAPTER 2, METHOD OF ANALYSIS 


2.1 Camera calibration 

The exact orientation and lens position of the camera 
were determined using the measured positions in the film 
image of the photomarkers and other surveyed objects, such 
as pressure gauges. This calibration procedure has been 
described in detail by Dewey et al (1975). The positions of 
all fiducial markers, as they appeared in a chosen calibra- 
tion frame, were transformed to a camera object plane in 
real space using the camera position, two camera orientation 
angles, and two "optical" parameters: a magnification factor 
and an axial rotation angle. The orientation angles and 
Optical parameters were calculated so that the two reference 
point images in the calibration frame, Pl and P2, transformed 
to coincide with their positions in the object plane calculated 
directly from the field survey data. Since an exactly surveyed 
camera position was not provided, the three camera position 
coordinates were also calculated, together with the two orienta- 
tion angles and the two optical parameters. The camera 
position "optimization" was accomplished by requiring that a 
third reference point image, P3, transformed to its object 


plane position as calculated from the survey data, and that 
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the distance between a pair of reference point images, P4 and 

P5, transformed to the corresponding distance in the object 

plane calculated from the survey data. The calibration 

procedure was iterative, and ended when the results of successive 
calculations differed by less than a pre-set tolerance. 

The field survey data used in the calibration procedure 
for Shot 10 are listed in Table 1. The calculated camera 
position (+ 0.1 feet) and orientation angles (+ 0.001 °) 
are listed in Table 2, together with the positions of all the 
calibration points in the object plane and their "shifts", 
that is, the differences between the transformed calibration 
points and their object plane positions calculated directly 
from the survey data. 

The object plane positions and shifts of the calibration 
points are also shown in Figure 3. Circles represent the 
transformed calibration points and squares represent positions 
calculated using the survey data. The larger circles indicate 
the three points that were made to coincide exactly, viz. 

Pl = Wl, P2 = W3, and P3 = 300Wl. The large half-circle 
indicates the final reference point (P5 = 300W2) used in the 
camera position optimization (P4 was set to equal P3). Cl and 
C2 denote the two charges. Frame centre is near the camera 
centering photomarker, VP2A. The photomarkers VP1A and VP1B 


show large shifts probably because they were not re-surveyed 
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after the pole on which they were mounted was lowered to replace 
the marker missing on Shot 8. No similar explanation could be 
found for the vertical shift of the foreground photomarker 
300W2. 

Consistency in calculating the optimum camera position 
from experiment to experiment was one of the criteria used 
in judging whether or not the optimum position was acceptable. 
The optimum camera positions used were the ones showing the 
least movement between experiments, on the assumption that 
the actual camera position was the same for each experiment. 
The calculated camera positions are shown below for the four 
experiments. Ground zero was given coordinates 2000 ft east, 
2000 ft north, and 2316.32 ft elevation. 


CALCULATED CAMERA POSITIONS 


East North Elevation 


2003.38 1385.09 2341.84 ft 
2002.76 2385.15 2341.66 
2002.62 1388.29 2341.72 


2003.25 1386.89 2340.64 


Range 0.76 3420 1.20 


The greater range in the north coordinates is probably due 


to the method of optimization used for this coordinate 


direction. The method was based on the apparent distance 


between two points in the object plane rather than on the 
apparent position of one point in the object plane. 

The position and orientation of a camera, once calculated, 
were assumed to remain constant throughout an entire film. 
The magnification factor and axial rotation angle were re-cal-~ 
culated for each frame in the data reduction part of the 
analysis described below, to account for frame-to-frame 


variations in the film projection optics. 


2.2 Data reduction 


The positions of the smoke puffs and two reference photo- 
markers as they appeared in the projected film image were 
digitized for each frame over a section of film that corresponded 
approximately to the duration of the positive phase of the 
blast waves. For the Shot 10 film, this time interval 
corresponded to frames 10 through 200 (frame 0 corresponded 
approximately to detonation time). Before each digitizing 
session, the projection system optics were aligned using a 
precision grid in the film gate of the projector. 

The two reference markers used for the reduction of Shot 
10 were P6 = VP3B and P7 = 300W2. The digitized data were 
transformed mathematically to the camera object plane using 
the image positions of these two reference markers in each 
frame and their positions in the calibration frame, to correct 


for frame-to-frame variations in image magnification and axial 


Lg 


rotation in the film projection optics and to correct for 
frame-to-frame effects of camera vibration. 

The camera object plane was defined to coincide as 
nearly as possible to the actual plane of smoke puffs. Smoke 
puff positions in the object plane were then transformed to 
the actual smoke puff plane, which was assumed to be a 
vertical plane running from ground zero at an angle of 6.7° 
south of west. The coordinate system in this plane was 
defined to have its x-axis horizontal, its y-axis vertical, 
and its origin at ground zero. Thus the entire grid of smoke 
puffs was positioned in real space for each frame. 

A time (after charge detonation) was assigned to each 
frame by computing the film speed at selected frames. Film 
speeds were computed using the 1 ms timing marks placed on 
the film at the time of its exposure. The reciprocal of a 
simple polynomial function fitted to the computed film speeds 
was then integrated to determine the successive frame times. 
Zero time was determined using the detonation zero timing 
mark placed on the film. This method of film timing has been 
described in detail by Dewey et al (1975). A complete set of 


timing data for the Shot 10 film is provided in Table 3. 


2.3 The output data plane 


The coordinate origin in the smoke puff plane was 
nominally ground zero. However, the coordinate origin was 


not located at the ground zero position which was listed with 
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the field survey data. The location used was instead a value 
corrected to account for the fact that the charges were not 
positioned in any experiment exactly above the surveyed ground 
zero. The smoke puff plane origin was defined to be that 
point which had the same elevation as the surveyed ground 
zero, but which lay directly below a point halfway along the 
line joining the charge centers. This corrected ground zero 
for Shot 10 was displaced horizontally by 1.2 ft from the 
surveyed ground zero, in a direction which was approximately 
30° south of west. 

Particle trajectory data were output in the smoke puff 
plane using an x-coordinate axis which was defined to run 
horizontally outwards from the corrected ground zero and 
a y-coordinate axis which ran vertically upwards from this 
same point, so that increasing x-coordinate values denoted 
increasing distance from an idealized vertical charge axis 
and increasing y-coordinate values denoted increasing distance 
above an idealized horizontal ground surface. All data in 
the smoke puff plane as they are presented in the figures of 
this report have x-coordinate values increasing to the right, 
so that the smoke puffs appear to run to the right of the 
charges rather than to the left as they appeared in the original 
photographic images. 

Figure 4 shows the positions of the 108 detonated smoke 


puffs comprising the grid for Shot 10 at a time prior to the 


eal 


detonation of the two charges. These positions are in the plane 
of the smoke puff grid and the charges, as described above. 

The smoke puff plane is not exactly parallel to the camera 

image and object planes (Figures 2 and 3), and various 
geometrical corrections were applied to make the transforma- 
tion between them. The puffs enclosed in parentheses were not 
visible in the earlier film frames, but were seen later when 
they were illuminated by the light of the fireball. Charge 
positions in the figures are plotted as if they were positioned 
exactly above the corrected ground zero origin. The data 


shown in Figure 4 have not been scaled. 


2.4 Scaling and ‘units 

The position-time histories of individual smoke puffs 
were extracted from the frame-by-frame positions of the smoke 
puff grid. This set of particle trajectories was then scaled 
to standard atmospheric conditions and charge weight to 
remove experiment-to-experiment differences due to variations 
in atmospheric conditions and charge weight. A change from 
British units to SI units was made at this stage of the 
analysis. 

Particle trajectory data were scaled by dividing all 
distances by Sachs scaling factor S = ¥ (WE )/ LP) and 


multiplying all times by the factor C/ (C58) , where C is 
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the ambient sound speed computed for Shot 10. Data used to 
compute C and S, and define the scaled event, are listed 


below with the computed values of C and S. 


Ambient temperature, 
Ambient pressure, 
Relative humidity, 
Computed vapour pressure, 


Ambient sound speed, 


Charge weight, 


Sachs scaling factor, 


Standard charge weight, 
Standard pressure, 
Standard temperature, 


Standard sound speed, 
(dry air) 


-5.94 °C 
94.38 kPa 
81.0 & 

O.S2 kPa 


328.003 m/s 


489.9 kg 


8.0718 


1.0 kg 
101.325 kPa 
ES 7'C 


340.292 m/s 


(223) 22) 


(23.689 PSt) 


(2.4 mm Hg) 


(1076 £t/s) 


(1080 lbs) 


(2.2 Ibs) 
(14.7 PST) 
(S59. OR) 


(L116 £t/s) 


The results presented in this report therefore apply to 


a scaled event which is the detonation of two 1 kg charges 
in a standard atmosphere. The scaled heights of burst for 


Shot 10 were 0.563 meters and 1.713 meters. 


2.5 Trajectory fitting 


Figure 5 shows the scaled particle trajectory data for 
Shot 10 in the smoke puff plane with positions measured 
horizontally and vertically from corrected ground zero. 
Approximately 9600 puff positions are represented. As 
represented, the raw trajectory data have not been edited, 
and a number of obvious film reading errors can be seen. 

The raw particle trajectory data were edited to remove 
obvious data prececeina erpars such as a single point widely 
displaced from its trajectory for one or two frames. The 
trajectories were then smoothed by fitting, for each puff in 
turn, simple polynomial expressions separately to both the x- 
and y-coordinate data, these being discrete functions of frame 
time. A least squares fitting technique was used. The adequacy 
of each fit was determined by examining on the same graphical 
output plots of both the raw trajectory data and the fitted 
curve. For Shot 10 this meant examining and adjusting 216 
such plots, at least two or three times each. 

For a given puff, the first step in fitting the raw 


trajectory data was to set the time of arrival of the shock 


] 
| 


front first hitting the puff. This was done subjectively, 


using the plotted raw data. The obvious constraint was 

used that this time for a single puff should be the same for 
the x- and the y-coordinates. Prior to this time of arrival, 
the measured positions of the puff were averaged to establish 


the initial puff position. 
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The raw trajectory data at times subsequent to the time 
of arrival were fitted with polynomial functions as described 
above. The fits were individually adjusted by varying the 
order of the polynomial and/or the relative weight given to 
each of the raw data points. The lowest order of polynomial 
which appeared to describe the particle trajectory was used. 
In all cases polynomials of fifth order or less were found to 
be adequate. Data points in obvious error were edited out by 
weighting them zero. For some puffs undergoing relatively more 
complex motions (i.e. those nearer to the charges), trajectories 
were fitted in two sections. If it appeared that the change 
in motion of the puff was the effect of a second shock wave, 
then it was not required that the first derivative of the 


fitting function (the velocity of the puff) be continuous. 


Otherwise it was required that this derivative at least appear 


to be continuous. In judging the motion of a particular puff, 

the motions of its near neighbours were also taken into account. 
Using the fitted functions, the positions of the smoke 

puffs were interpolated at a series of discrete times, which 

were not necessarily the frame times. The first derivatives 

of the fitted functions were also calculated at the interpolated 


times for use in later calculations of particle velocity. 


2.6 Regionalization and shock strength calculations 

Knowing the time of arrival of a shock front at each 
puff, and the position of each puff, calculations were made 
of the shock trajectory, the shock velocity, and various 
peak parameters associated with shock front strength. These 
shock trajectories and associated peak parameters were then 
compared to the corresponding results calculated as part of 
the refractive image analysis already reported (Dewey et al, 
LOTS) 

To associate a shock radius with a given puff position, 
the smoke puff plane was divided into five regions. All 
puffs in a given region were assumed to have been first hit 
by the same shock front. The region boundaries were defined 
using triple point path measurements obtained from the refractive 
image analysis mentioned above. 

The regions which were defined for Shot 10 are shown in 
Figure 6. In regions 1 and 2 the puffs were assumed to have 
been first struck by primary shock fronts, and the shock radii 
were measured from the appropriate charge center. Puffs in 
regions 3, 4 and 5 were assumed to have been first hit by a 
Mach stem: along the ground, and below and above the interaction 
plane, respectively. The interaction surface, along which the 
upper and lower primary shock fronts interact without energy 
loss, was assumed to be the plane of symmetry separating the 


charge centers. 
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In the two Mach stem regions adjacent to the interaction 
plane, radius values were meaSured from the point halfway 
between the charges, and in the Mach stem region adjacent to 
the ground surface they were measured from corrected ground 
zero. This was done in accordance with the assumptions that 
the Mach stem shock fronts were spherical and that the sphere 
centers remained fixed with time at the points where the vertical 
coordinate axis intersected the interaction and ground planes. 
Shock strength results computed in this way for Shot 10 represented 
an improvement over similar results computed using the assump- 
tion that the Mach stem fronts were cylindrically shaped and 
centered about the vertical coordinate axis. This question 
of Mach stem shape is discussed further in a later section of 
this report. 

In each of the five regions, the shock trajectory data 
obtained from the first movement of the smoke puffs were 


fitted to a function of the form 


Ete) = A+ BE +E tog(i + =) ; 


where r is the shock radius, t is the time after detonation, 
and A, B, and C are the fitted coefficients. The shock 
velocities were calculated by differentiating this function. 


and peak 


The peak particle velocity, wy peak density, Der 


hydrostatic overpressure, Por in each of the five regions were 
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calculated as functions of shock radius, using extensions of 


the Rankine-Hugoniot equation, viz: 


Ys. 2 M?-1 
2 

Cc. yt+l M 

Ps _ _(y+1)M?_ 


Ds (y-1)M*+2 


2Y (y2- 
P yt+l (M"-1) 


where Ca» B, and ¥ are the ambient values of sound speed, 
density and pressure, respectively; y is the ratio of specific 
heats and M the Mach speed of the shock front. Peak particle 
velocity will be used in this report as a measure of shock 
strength because other independent calculations of this para- 
meter were also made. Peak values of density and pressure 


were calculated using the above method only. 


2.7 Particle velocity calculations 


Assuming the smoke puffs to be perfect tracers of the 
gas flow within the blast waves, the gas particle velocities 
throughout the field of smoke puffs were calculated for a 
sequence of times. This was done by calculating the derivative 


of the polynomials used to interpolate the fitted trajectories. 


The derivatives of these polynomials were also calculated at 


the times-of-arrival of the shock at each puff to provide a 
second measure of V,, the peak particle velocity immediately 
behind the shock. 

The peak particle velocities ¥.. were also calculated by 
a third method. At specific times the particle velocities, 
V, and the radial positions, r, of all the puffs within each 


of the five regions were fitted to a function of the form 


that is, a straight line with positive slope. (The data 

scatter did not justify the fitting of a more flexible function.) 
This fitted function, at each time and in each region, was 
extrapolated to the appropriate shock radius at that time 


to determine peak particle velocity. 


- 


2.8 Density and hydrostatic overpressure calculations 


The particle trajectory data were also used to compute 
time varying air densities throughout the field of the smoke 
puff grid. The calculated densities were in turn used to 
compute hydrostatic overpressures in the same field. 

The smoke puff grid of 9 columns of 12 puffs each was 
treated as a grid of 88 cells, each cell being defined by 


4 adjacent smoke puffs. Each such quadrilateral cell was 


29 


viewed as representing the cross-section of an elemental 

volume defined by rotating the area of the cell about the axis 
of cylindrical symmetry rising vertically through the two 
charges. Conservation of mass within this volume demands that 
the product RAD is constant at all times, where R is the radial 
coordinate of the centroid of the cell defining the volume, 

A is the area of the cell, and D the average density of air 

in the volume. In particular, the product RAD must always 
equal (RAD) ,, its ambient value computed before the arrival 

of any shock fronts at the cell. This equality gives the 


relative density D/D, at any time, as 


Pressure values were calculated using the fact that p/p’ 
remains constant as long as the entropy of the gas in a cell 
remains constant, where P is absolute pressure at a point, D 
density, and y the ratio of specific heats of air. In 
particular, the, value of P/D’ equals the value py mat 
measured immediately behind the leading shock front, and remains 


constant until the entropy is changed by the arrival of a 


second shock. The relative overpressure at a point can be 


computed using 


The values ee and D/D, were computed at each puff's 
position immediately after it was first struck by a shock 
front, taking r. and Ds from the shock velocity data computed 


according to the method described earlier (section 2.6). The 


values computed for PSF and D./D, at the cell corners were 


averaged. The overpressure P/F. - 1, subsequently calculated, 
therefore also represents an average value over a given cell. 
The arrival of shock fronts subsequent to the first could not 

be detected accurately, so that no attempt was made to determine 
exactly those times after which the calculated overpressures 


in given cells became invalid. 


2.9 Summary 


Previous sections describe the method of particle trajectory 
analysis by which data of four basic types were obtained: times 
of shock front arrival, particle velocities, densities and 
hydrostatic overpressures. To obtain a measure of shock 
strength the particle trajectory data were analyzed in separate 
groups according to the shock front which first arrived at 
the smoke puffs in each group. It was assumed that the 
shock fronts had a specific regular shape. In this way the 
smoke puffs at the observed times of shock front arrival were 
assigned radial position coordinates to supplement their 


measured Cartesian coordinates. Radial shock velocities 


throughout the smoke puff grid were computed and peak 
particle velocities were derived from the shock velocities. 
These peak velocity values were used as a measure of shock 
front strength and were compared with values obtained using 
smoke puff data in two other ways, and with values obtained 
previously using refractive image analysis. ‘he particle 
trajectory analysis depended on the refractive image 
analysis only to define the five regions in the smoke puff 
plane. ‘These regions were used only in the analysis of the 
shock fronts and were of no particular significance in the 
subsequent analysis of the smoke puff trajectories behind the 
shocks. 


The basic data obtained from the smoke puff trajectories 


after shock front arrival were derived in Lagrangian coordinates, 


that is, for a grid of points which moved with the air in the 
blast wave and did not remain fixed with respect to the 
ground surface. The method used to transform these data to 
Eulerian coordinates, that is, fixed in space, and which 

also permitted the calculation of dynamic pressure fields, 


is the subject of the next chapter. 
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CHAPTER 3. 


SURFACE REPRESENTATION 


3.1 Introduction 


Measurement of the smoke puff positions permitted the 
determination of several physical properties of the blast 
waves in a vertical x-y plane passing through the corrected 
ground zero. Letting z denote values of a particular physical 
property at a specific time, the results of the particle 
trajectory analysis may be viewed as sets of discrete points 
on three-dimensional surfaces, f(x,y,z) = 0. If such surface 
representations are made explicit, they may be used to inter- 
polate z-coordinate values, i.e. values of the physical 
property, at x-y positions other than the actual smoke puff 
positions, such as gauge locations. Such surface representa- 
tions also may be used to display graphically an entire set of 
physical properties in the form of contour lines joining points 
of equal z-value. The following sections describe how the 


various surfaces were fitted. 


3.2 Preliminary interpolations 


Source data (z-coordinates) for surface fitting were the 
times-of-arrival of the first shock front at the individual 
smoke puffs, the measured particle velocities at puff positions 
at a sequence of times 


(velocity fields), and average densities 


at cell centroids at the same sequence of times (density fields). 
Values of hydrostatic overpressure at the cell centroids 

were also fitted, these values being derived from the density 
data, as described previously. The x-y coordinates were 

the initial puff positions, the puff positions at a time, 

t , and the cell centroid positions at a time, t . 

The first step in the surface fitting procedure was 
the completion of the data grid. This meant generating x-y 
coordinates to represent "missing" puffs and cell centroids, 
and interpolating z-coordinates at these new positions. This 
was done in a very conservative manner, using linear inter- 
polation in a domain bounded by the completely convex set 
of straight line sections which connected the outermost of the 
Originally existing puffs. No extrapolation beyond this 
domain was attempted. 

The positions of missing puffs were assigned by inter- 
polating x- and y-coordinates in turn, between those of the 
adjacent puffs in each missing puff's row and column. The 
z-data assigned to the new positions were then interpolated 
along the row and along the column, and the two results were 
averaged. New z-coordinates were assigned only to new puff 
locations that fell inside the original domain boundaries. 

A new puff position on the domain boundary, that is a puff 
with only one nearest neighbour in its row or column, was 


given a position coordinate and a z-coordinate equal to 
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those of its single neighbour in the appropriate row or 
column calculation. 

In Shot 10 there were no missing puffs in the original 
grid and it was not necessary therefore to apply the above 
technique when fitting the time of arrival surface. In 
all other cases, missing puff data were generated only to 


expedite subsequent interpolations. 


Fas Interpolation onto an Eulerian grid 


The second step in the surface fitting procedure was 
the definition of a rectangular grid and the subsequent 
interpolation of a z-coordinate value at each new grid 
point inside the original domain. In all cases the new 
grid was a regular square grid with a constant mesh size 
of 0.1m, which was approximately one-half the mesh size 
of the original smoke puff grid after scaling. The new 
grid was in the same plane as the smoke puffs but, unlike 
the smoke puff grid, it remained stationary relative to 
the coordinate origin. 

Interpolation at the new grid points was done by fitting 
a polyhedral surface to the original, completed smoke 
puff grid. ‘This was done by fitting plane surfaces to 


triplets of data at the corners of the original smoke puff 


cells. For example, when a new grid point fell inside an 


original cell of smoke puffs (1,2,3, and 4) as shown: 


new grid line 


a z-coordinate value was assigned to the new grid position "6" 
by flitting planes to the puff data at 1, 2, and 3, and at 2, 3, 
and 4, and interpolating at the new grid point using each 
plane in turn, and averaging the two results. If, at the 

edge of the domain, a cell was defined by only three puffs, 
only one interpolation plane was used. When two or more cell 
corners were beyond the domain, no data were interpolated 

to the new grid. In this way the original domain boundaries 
were Shrunk, and made concave, on occasion. 

The surface represented by the resulting grid of inter- 
polated data was therefore an array of connected triangular 
plane elements,with four such elements for each complete 
smoke puff cell. The connected array was everywhere continuous 
in z, albeit discontinuous in the derivatives of z across 
the connections, and it passed through all of the original 


smoke puff data. There was therefore no smoothing of the 
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original data, and no chance of excessive values of z between 


original data due to uncontrolled fitting functions. 


3.4 Dynamic pressure calculations 

Surfaces were fitted to the times of shock arrival at 
each smoke puff and to the particle velocities, densities and 
hydrostatic overpressures at a sequence of times. The grids 
were Eulerian in that their oosition in the x-y plane did 
not change with time. The z-coordinates were, specifically, 
the velocity vector magnitudes expressed as Mach numbers, 
density ratios relative to ambient density, and overpressure 
ratios relative to ambient pressure. 

Since the same interpolated grids were used for velocity 
and density, the dynamic pressure at each grid point also 


could be computed using 


be | i Ce / | | é 
dynamic pressure ratio = 2eiv) a =) @xak 
P ZN | Cc 
fo) fe) fe) 
where y = 1.4 ,; (D/D,) is the measured relative density, 
(|v[/c,) the measured particle velocity Mach number, e 


the standard pressure, and o. the standard sound speed. 
Dynamic pressures could not be calculated prior to this point 
in the analysis because the original, uninterpolated values 
of relative density and particle velocity were not measured 


at the same set of points in the x-y plane. 
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3.5 Contours and time histories 

The fitted surfaces described above were used only for 
interpolation; the data were not smoothed. Interpolation 
onto an Eulerian grid made possible the calculation of 
dynamic pressure fields. Further interpolation to points 
inside the Eulerian grid made it possible to compute and 
plot contour lines over the various data fields, and also 
time history curves for data at specific fixed locations within 
those fields. 

The contours calculated were those of equal time of 
shock arrival (isochrones), particle velocity (isotachs), 
density (isopyknics), hydrostatic overpressure (static 
isobars), and dynamic pressure (dynamic isobars). The 
isochrones represent the shock front shapes at different 
times. The contouring technique therefore provided a 
complete mapping of all the physical properties of the blast 
wave within the smoke puff domain without reference to other 
measurements. The time histories of these physical properties 
at several fixed locations were also calculated. Some of 
these were calculated at pressure gauge positions in unscaled 
units, so they may be compared to actual gauge measurements. 

Contours of specific z~values were calculated using | 
linear interpolation, searching between pairs of grid points 
along the Eulerian-Cartesian grid lines and along the mesh 


diagonals for possible contour line intersection points, and 
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joining the points which were found. Time history curves 


were calculated by interpolating z-values at a fixed location 
in a time-series of data fields which had been interpolated 
onto the same Eulerian grid. The interpolations to points 
inside this common grid were of the polyhedral surface 

type describe previously. Because all fitted surfaces were 
continuous in z and because the mesh size could be made 
relatively small, satisfactory contours and time histories 
were obtained. Because the fitted surfaces were smooth 
(comparitively speaking) and, in the limit of an infinitely 
small mesh, they would pass through all the control points in 
the measured data fields, the contours and time histories 
were smooth and represented the measured data faithfully in 


all cases. 


CHAPTER 4. SHOT 10 RESULTS 


4.1 Times of shock front arrival 

The measured initial puff positions, the times of first 
shock arrival, and the peak particle velocities obtained by 
differentiating the functions fitted to the particle 
trajectories are presented in Table 4. Puff position is 
given relative to corrected ground zero as origin with 
horizontal and vertical axes. Puff position and the time 
of arrival of the first shock are given both as observed and 
scaled. Particle velocities listed are derivatives of the 
fitted puff trajectories at the times of shock arrival, and 
are expressed in Mach units. Expressed this way, the particle 
velocities are the same scaled as unscaled. Also listed are 
the initial radial puff positions (scaled) and region codes. 

Shock front data determined from the first movement of 
the smoke puffs, i.e. calculated from the time-of-arrival 
data in Table 4, are listed in Tables 5.1 - 5.5. Each 
table corresponds to one of the 5 regions used. Listed are 
the observed and fitted unscaled shock trajectory data, the 
scaled fitted shock trajectory data, and the computed shock 


velocities and peak parameters associated with shock front 
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strength: peak hydrostatic overpressures in atmospheres and 
in kilopascals, peak particle velocities in Mach units, and 
peak density ratios. Given as ratios, these peak parameters 
are the same scaled as unscaled. Pressure given in kilo- 
pascals in the tables refers to the unscaled (observed) 

case only. 

The shock front radius versus time data derived using 
particle trajectoryanalysis (PTA) are also shown in Figures 
7.1 - 7.3 for the two primary fronts, the two Mach stems 
at the interaction plane, and the ground Mach stem, 
respectively. They are compared to corresponding data 
derived from refractive image analysis (RIA) reported by 
Dewey et al (1975). The refractive image analysis results 
were obtained using photogrammetry against a striped canvas 
backdrop and they describe the shock as it travelled ina 
direction almost diametrically opposite to the direction 


of the smoke puff grid. 


4.2 Shock strengths 
Peak particle velocities calculated from shock front 
velocities are shown in Figures 8.1 - 8.3 for the primary 
fronts, interaction Mach stems, and the ground Mach stem. 
This method of determining peak particle velocities has 


been labelled method 1, and the data plotted correspond to 
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those listed in Tables 5.1 - 5.5. The results in the figures 
are compared with those previously obtained using refractive 
image analysis (RIA). In the case of the primary shock fronts, 
results are also compared to those of Brode (1957) for TNT. 

Peak particle velocities calculated as the derivatives 
of the fitted particle trajectory functions at the time of 
shock arrival are shown in Figures 9.1 - 9.3, labelled method 
2. Corresponding numerical data are listed in Table 4. 

Results are compared to Brode in the primary front case, and 
with the results of method 1 in the Mach stem cases. 

Peak particle velocities calculated by fitting particle 
velocity as a function of radial position in a region at a 
fixed time, and extrapolating to the shock radius at that 
fixed time (method 3), are shown in Figures 10.1 - 10.3. 

The corresponding numerical data are listed in Tables 6.1 - 
6.5. Peak particle velocities are compared to Brode in the 
primary front case, and to method 1 results in the Mach stem 
cases. 

Some of the data from which the method 3 results were 
extrapolated are shown in Figure 11.1. Shown are the magnitudes 
of the particle velocity vectors throughout the blast wave in 
the regions of the Mach stem regions above and below the 
interaction plane. Results are plotted versus a radial 


coordinate computed in the same manner as for the time-of- 


arrival results. Also shown are the straight line fits to 


the data, extrapolated to the shock front radius at each of 
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the times shown. Figure 11.2 shows similar results for the 


primary shock front regions, compared to the results of 


Brode for TNT. 


4.3 Particle velocity fields 


The calculated particle velocities in the plane of the 
smoke puffs are shown as vectors in Figures 12.1 through 
12.23 for times between 0.4 and 1.0 ms in 0.1 ms increments, 
and for times between 1.0 and 6.5 ms in 0.5 ms increments. 
The particle velocity field is shown at several other inter- 
mediate times to illustrate certain columns of smoke puffs 
being struck by the advancing shock fronts. All times and 
positions are scaled to a 1 kg charge in a standard 
atmosphere. Position coordinate axes extend horizontally and 
vertically from the corrected ground zero. The charge 
positions are also shown. The particle velocity vectors 
represent the derivatives of the smoothed particle trajectories, 
and their magnitudes may be judged using the standard vector 
shown on each figure. All velocities are measured in Mach 
units, relative to the standard sound speed. Puffs not yet 
struck by a shock wave are represented by small circles 
(zero velocity). 

Numerical data corresponding to Figures 12.1 - 12.23 are 
listed in Tables 7.1 through 7.16, along with scaled radial 


positions of the puffs, and region codes as defined in Figure 6. 
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Conversion factors are given at the foot of each table, which 
may be used to convert the scaled data in the tables and 


figures back to their original unscaled values. 


4.4 Density and hydrostatic overpressure fields 


Calculated average relative densities throughout the smoke 


puff plane are depicted graphically in Figures 13.1 - 13.23, 


for times between 0.4 and 1.0 ms in 0.1 ms increments, and 


between 1.0 and 6.5 ms in 0.5 ms increments. 


The average 


density field is shown at several other intermediate times to 


illustrate the transition of certain cell columns through 


their periods of peak density. All time values are scaled. 


Cell positions are scaled and are given relative to the corrected 


ground zero as origin with horizontal and vertical axes. Charge 


positions are also shown. The calculated densities may be 


judged using the density shading scale shown on each figure. 


Density is given as a ratio, relative to ambient density. 
Cells not yet struck by a shock wave and cells in which the 
density has dropped to a value less than ambient density are 
shown blank. 


Corresponding numerical data are listed in Tables 8.1 - 


8.8 along with radial cell positions computed according to the 


regions defined previously. 


Numerical data describing the 
fields of hydrostatic overpressure are Similarly listed in 


Tables 9.1 - 9.8. The pressure results for a given cell were 


obtained by multiplying the density results for that cell by 
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a factor determined by the strength of the shock which first 


traversed the cell, but which did not vary with time. 


4.5 Times-of-arrival surface 

Figure 14 shows a perspective view of the surface fitted 
to the original smoke puff positions and the observed times 
of first shock front arrival, i.e. to the data listed in 
Table 4. The grid mesh size is 0.1 by 0.1 meters (scaled), 
about 2.5 feet square (unscaled), or about 1/2 that of the 
original smoke puff grid. The charge positions are indicated 
on the vertical distance axis. 

The times-of-arrival surface is smooth enough to permit 
contouring and the contours in this case (isochrones) represent 
shock front shapes at different times. These contours are 
plotted in Figure 15. The times-of-arrival surface was not 
smooth enough to permit the calculation of gradient vectors which 
could be used to compute shock velocity vectors and shock 
strengths over the new grid. 

Two attempts were made to obtain contours of shock strength. 
In the first, the times-of-arrival surface was smoothed by 
least-squares fitting low-order, one-dimensional polynomial 
functions to the time-of-arrival data along each grid row and 
column separately, and computing the derivatives of the fitted 


functions to obtain the associated components of the surface 
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gradient vectors. Shock velocity vectors were obtained from 
the time-of-arrival gradients, and from these peak particle 
velocities were computed. The peak particle velocity (shock 
strength) surface is shown in Figure 16. The contours of 

this surface (not shown) did not exhibit any discontinuities 
across the boundaries of the shock front regions, as they would 
if surfaces were fitted to the times of arrival in each region 
separately. 

The results of the second method used to compute shock 
strength contours are shown in Figure 17. These were obtained 
by interpolating an average shock radius at each value of 
peak particle velocity shown, for each shock front region in 
turn, using a composite of the various peak particle velocity 
versus radius curves shown in Figures 8 through 10. Arcs of 
circles with these radii, centered on the appropriate points 
along the vertical charge axis, were then drawn in the regions 
to represent shock strength contours. These peak value contours 
are discontinuous across triple point locii and other region 
boundaries. As.a result, some horizontal lines are crossed 
twice by the same contour or, in other words, identical shock 
strengths can be found at two locations the same vertical 
distance from a reflecting surface, but at different radial 
distances from the vertical charge axis. Examples are illus- 


trated in the figure. 
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4.6 Field surface contours 

Contours of equal particle velocity, density, hydrostatic 
overpressure, and dynamic pressure in the blast waves were 
determined for a series cf times, using surfaces fitted to the 
various measured data fields at those times. Sample results 
are shown in Figures 18 through 21 at scaled times of 2.5 ms 
and 4.0 ms. The shock fronts shown in these figures are obtained 
from the time-of-arrival surface (as were those in Figure 15). 
Field contours such as those shown can be drawn for any scaled 
time between 0.5 ms and 6.5 ms. It should be re-stated that 
all of these results were obtained from the photography of the 
smoke puffs only and do not rely on the results obtained using 


the refractive image analysis (Dewey et al, 1975). 


4.7 Time histories 

By mapping the physical properties of the blast waves at 
short time intervals it was possible to determine the time 
histories of these properties at any selected fixed position 
within the smoke puff grid. This was done at 12 fixed locations, 
three in the two primary regions and three in each of the three 
Mach stem regions, as shown in Figure 22. At each distance from 
the axis of the charges in the Mach stem regions, each of the 
time history stations is the same distance from either the 


interaction plane or the ground plane. Particle velocity time 


histories could be interpolated closest to the ground level 


because these were measured at puff locations, whereas the 
density and pressure data were measured at cell centers. 

Time histories of particle velocity, density, hydrostatic 
and dynamic overpressure at these locations are given in 
Figures 23 to 26. 

Also plotted with the time histories are the interpolated 
values of the time of arrival of the first shock front at the 
stations. The height of this time-of-arrival line represents 
a peak value derived from the shock velocity analysis. 

Time histories for hydrostatic and dynamic pressure are 
also plotted in Figure 27.1 to 27.7 for stations at the 
nominal positions of field-mounted pressure gauges on the 
"60 foot gun barrel". The gauges on this gun barrel were 
mounted at nominal elevations of 10, 15, 20, 27, 30, 33, and 
40 feet. The time histories at these locations are given in 
unscaled units in order to facilitate comparisons with the 
gauge measurements. 

The dynamic pressures plotted in figures 26 and 27 are 
Maximum values, computed using both the x and y components of 
particle velocity. Similar plots were made of the horizontal 
components of dynamic pressure, but the differences were not 
significant since the y components of particle velocity at these 
locations were small. Other locations could have been chosen 


at which the y components would not have been insignificant. 
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CHAPTER 5. DISCUSSION 


5-1 Choice of Dipole West Shot 10 


The results presented in this report were obtained from 
an original attempt to analyze an array of discrete two- 
dimensional particle trajectories within a blast wave. 
(During the preparation of this report Suffield Special Publica- 
tion No. 71, by John Anderson was published, which also 
describes a blast wave two-dimensional particle trajectory 
analysis using somewhat different techniques to the ones 
presented here.) Previous analyses of blast wave particle 
trajectories were limited to those cases in which it could 
be assumed that the motion of the particles was rectilinear 
(Dewey, 1964; 1971). 

Particle trajectory photogrammetrical measurements were 
made on Dipole West Shots 8, 9, 10 and 11, and it was decided 
to develop the two-dimensional analytical techniques using 
the measurements from Shot 10 because in this experiment a 
smaller number of smoke puffs had been used, and these puffs 
had formed with a high degree of success so as to produce an 


almost always complete rectangular array of particle tracers. 


5.2 Particle trajectory analysis technique 
The simplest possible techniques were used to analyze 


the particle trajectories obtained for Shot 10. The trajectory 
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of each puff was described by least squares fitting a low order 
polynomial to the time variation of both the x and y coordinates 
of the puff. The two polynomials obtained in this way could be 
differentiated at any point to give the particle velocity 
components, and interpolated to give the particle position at 
any time. The gas density within the blast wave was determined 
by using the positions of four adjacent puffs to define a 
quadrilateral cell. The resulting density is thus an estimate 
of the average density within the cell, the position of which 
is described by its centroid. 

This analysis system considers each smoke puff in isolation, 
and although some smoothing of the data is achieved by the 
least squares fitting, no account is taken of the positions of 
adjacent puffs. Subsequently linear interpolation is used 
between the smoke puff positions so as to obtain a complete 
mapping of the blast wave properties throughout the domain of 


the smoke puff array. 


Fes Reliability and accuracy of the measurements 


An attempt has been made in this report to evaluate the 
reliability of the particle velocity and density measurements 
obtained from the analysis of the time resolved smoke puff 
trajectories. This is not easy to do since there are no 
similar measurements with which the results can be compared. 


However, the particle trajectory results have been compared 
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with the peak shock front values obtained from the refractive 


image shock velocity measurements (Dewey et al, 1975). Also, 
hydrostatic and dynamic pressure results obtained from the 


particle trajectory analysis are presented for positions at 
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which pressure gauges were mounted so that a direct comparison 
of the pressure-time histories will be possible. 

Comparisons with the peak values at the shock front are 
not completely valid because the particle trajectory analysis 
method does not give results at the shock itself but only in 
the blast flow behind the shock. This necessitates extrapolat- 
ing the particle trajectory data to the expected position and 
time of arrival of the shock. Extrapolation is rarely reliable 
and in the technique described as method 3, extrapolation was 
applied to data gathered from a large region of the blast wave 
and which thus showed significant scatter. 

In order to calculate shock velocities assumptions were 
made about the shapes of the shock fronts. In the case of 
the primary shocks the assumption of sphericity appears to be 
valid. The shape of the Mach stem shock is still an open 
question. In the refractive image analysis (Dewey et al, 1975) 
the Mach stem shock positions were measured as close as 
possible to the reflecting surfaces and the assumption of 
cylindrical symmetry appears to have been valid. However, 


in the analysis of the smoke puff trajectories, measurements 
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could not be made so very close to these reflecting surfaces. 
It was clearly seen that the Mach stems were not cylindrical 
and as a first approximation it was assumed that these shocks 
were spherical with the centre at ground zero, or at the 
equivalent point on the ideal reflecting surface. This 
assumption gave results which were internally consistent, 
but it is clear that an investigation of the shape of the 
Mach stem shocks produced by a spherical shock reflecting 
from a plane surface, is a matter of urgency. 

Three methods of calculating shock strength are described. 
The first uses the time of arrival of the shock front at the 
smoke puffs within a given region and analyzes these data in 
the same way as for the shock front refractive image photo- 
grammetry described by Dewey et al, 1975. These results 
(PTA) are compared with the refractive image (RIA) results in 
Figures 8.1 to 8.3. There seem to be some significant 
differences between the results of the two measurement 
techniques: in particular the lowering of the Mach stem shock 
strength over the rough ground compared with the ideal reflect- 
ing surface is not so apparent when measured using the particle 
trajectory technique. This is probably due to the fact that 
the refractive image measurements were made within 0.3 m 
of the rough ground surface whereas the particle trajectory 


time-of-arrival measurements were made throughout the Mach 


stem region up to a distance of 4 m from the ground. Also, 


the two sets of results describe the shock fronts travelling 
in opposite directions from the charges. 
The second method of calculating peak particle velocities 


(shock strengths) used the derivatives of the two polynomials 


ee ee 


fitted to each particle trajectory, evaluated at the estimated 
time of shock arrival. Differentiation of a fitted polynomial 
at the extreme end of the observed data range is of questionable 
validity, and as might be expected the results plotted in 
Figures 9.1 to 9.3 show a considerable amount of scatter. 

This scatter is, however, uniformly distributed about the 
expected curve. 

The peak particle velocities were also calculated by 
fitting all particle velocities in a region as a function of 
radial position in the region at a fixed time, and extrapolating 
to the shock radius at that time. Because of the scatter in 
the data obtained throughout the region, fitting the results 
to a curve of higher than first order did not seem justified. 
The results of these extrapolations are shown in Figures 10.1 
bo Oa. 

Although the extrapolation of the particle trajectory data 
beyond its domain of validity has produced a large degree of 
scatter, the peak particle velocity results are generally 


consistent with the results calculated using the method of 


refractive image analysis. 
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An attempt has been made to estimate the inherent accuracy 
and reliability of the particle trajectory measurements. A 
smoke puff typically has a diameter of about 0.3 m, and it is 
estimated that the centre of such a puff is identified with 
an accuracy of approximately +0.1m. This leads to an estimated 
maximum error in particle velocity of +0.008 m/ms (+Mach 0.035), 
and an average density ratio of +0.35. The mean peak particle 
velocity and density ratio measured in this experiment were 
approximately Mach 1 and 2.5 respectively, and so it is estimated 
that the accuracy of the reported results is approximately 5% 
for particle velocity and 15% for density ratio. 

Dynamic pressure ratios were computed from the velocity 
and density results and are estimated to be accurate to within 
about 20%. Hydrostatic pressure ratios were computed using ] 
the density results and the peak overpressures computed from 


the shock velocities. Because the shock velocity data depended 


| 
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on region definitions and, in the case of Mach stems, on assumed 
shock front shapes, the accuracy of the hydrostatic overpressures 
computed behind: the shock fronts is the most difficult to 


judge. Hydrostatic pressure ratios computed using cylindrical 


and spherical Mach stem shapes differed by as much as 10%. 
Combining this maximum figure with the estimated uncertainty 
given above for densities, it is estimated that the inaccuracy 


of the hydrostatic overpressures reported may be as high as 25%. 
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5.4 Surface mapping and contouring 


The initial analysis of the smoke puff positions measured 
from the films provided a measure of the particle velocity of 
each smoke puff and the average density in each cell of four 
adjacent puffs. Chapter 3 described the use of linear and 
planar interpolation between the discrete values of particle 
velocity and density so as to provide a complete mapping of 
the parameters in the x-y plane, at any specified time. From 
such mappings, various contours were drawn to connect points 
having the same parameter value. In addition, from sequences 
of such mappings it was possible to determine the time 
histories of the parameters at several fixed positions within 
the smoke puff array. 

A considerable amount of information is obtained from the 
analysis of the particle trajectories, and it seems clear that 
contour mappings and time histories best display this information. 
The contours can be compared easily with similar results from 
computer code predictions of the blast properties, and the 
time histories of hydrostatic and dynamic pressure can be 
readily compared with gauge measurements. To facilitate future 
comparisons the contours can be mapped at any specified time 
during the passage of the blast wave, and the time histories 
can be calculated for any specified position within the smoke 


puff array. Such outputs will be provided on request. 
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5/5 Dynamic pressure corrections 


The dynamic pressure presented in this report is defined 
as one half the local gas density times the square of the 
maximum velocity component. In comparing these results with 
total head gauge measurements, compressibility corrections 


should be applied in the higher velocity regions. 
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